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The density of accessible levels at low spin in the ("^136,0:7) reaction has been extracted for 
the ^''^•I'^^Dy and i^i'i'^^y^ nuclei. The entropy of the even-odd and even-even nuclei has been 
deduced as a function of excitation energy, and found to reach a maximum of 15 fcs before neutron 
evaporation. The entropy of one quasi-particle outside an even-even core is found to be 1.70(15) 
fcs. This quasi-particle picture of hot nuclei is well accounted for within a simple pairing model. 
The onset of two, four and six quasi-particle excitations in the ^''^Dy and ^^^Yb nuclei is discussed 
and compared to theory. The number of quasi-particles excited per excitation energy is a measure 
for the ratio of the level energy spacing and the pairing strength. 

PACS number(s): 21. 10. Ma, 24. 10. Pa, 25. 55. Hp, 27.70. +q 



I. INTRODUCTION 

Already in 1936 H. Bethe introduced the Fermi gas to 
describe nuclear properties at high temperature Q . This 
simple independent particle model has been modified by 
including residual interactions between the nucleons. In 
the low excitation region long-range pair correlations play 
an important role and are roughly described within the 
so-called back-shifted Fermi gas model Q. 

There is evidence for the existence of paired nucleons 
(Cooper pairs) at low temperature. In high spin physics, 
the backbending phenomenon is a beautiful manifesta- 
tion of the breaking of pairs. The mechanism is induced 
by Coriolis forces tending to align single particle angular 
momenta along the nuclear rotational axis . Theoret- 
ical models also predict reduction in the pair correlations 
at higher temperatures |^^. 

The breaking of pairs is difficult to observe as a func- 
tion of intrinsic excitation energy. Recent theoretical 0] 
and experimental |^,^ works indicate that the process 
of breaking pairs takes place over several MeV of exci- 
tation energy. The corresponding critical temperature 
is measured to be Tc ~ 0.5 MeV/fcs ||l^], where fc^ is 
Boltzmann's constant. 

The aim here is to extract the entropy of the i6i,i62£)y 
and i7i,i72Yb isotopes, and deduce the number of ex- 
cited quasi-particles as function of excitation energy. In 
Sect. 2 we describe the experimental techniques and an- 
alyzing tools. Sect. 3 presents results for the entropy us- 
ing a simple pairing Hamiltonian with an even and odd 
number of fermions distributed over L single-particle lev- 
els with double degeneracy. This is a model, which for 
small numbers of fermions, typically less than 20, can 
be solved numerically yielding all possible eigenstates. It 
results in the full level density and can in turn be used 
to extract thermodynamical quantities. Since we expect 
pairing correlations to be important in nuclei, such a sim- 
ple model should mimic to a certain extent the entropy 



extracted from the experimental level density. In Sect. 4 
we present the experimental findings and relate them to 
the simple pairing model of Sect. 3. Concluding remarks 
are given in Sect. 5. 



II. EXPERIMENTAL METHODS 

The Oslo cyclotron group has developed a method to 
extract nuclear level densities at low spin from measured 
7-ray spectra l^-p^j. The main advantage of utilizing 
7-rays as a probe for level density is that the nuclear 
system is likely thermalized prior to the 7-emission. In 
addition, the method allows for the simultaneous extrac- 
tion of level density and 7-strength function over a wide 
energy region. 

The experiments were carried out with 45 MeV 
■^He-projectiles accelerated by the MC-35 cyclotron at 
the University of Oslo. The experimental data were 
recorded with the CACTUS multidetector array us- 
ing the (^He,a7) reaction on lea.ies^y i72,i73Yb self- 
supporting targets. The charged ejectiles were detected 
with eight particle telescopes placed at an angle of 45° 
relative to the beam direction. Each telescope comprises 
one Si AE front and one Si (Li) E back detector with 
thicknesses of 140 and 3000 /im, respectively. 

An array of 28 Nal 7-ray detectors with a total effi- 
ciency of ~15% surrounded the target and particle detec- 
tors. In addition, two Ge detectors were used to monitor 
the spin distribution and selectivity of the reactions. 

From the reaction kinematics the measured a-particle 
energy can be transformed to excitation energy E. Thus, 
each coincident 7-ray can be assigned to a 7-cascade orig- 
inating from a specific excitation energy. The data are 
sorted into a matrix of {E, E^) energy pairs. At each ex- 
citation energy E the Nal 7-ray spectra are unfolded Q , 
and this matrix is used to extract the first-generation (or 
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primary) 7-ray matrix with the well-estabHshed subtrac- 
tion technique of Ref . . 

The resuhing matrix P{E,E^), which describes the 
primary 7-spectra obtained at initial excitation energy 
E, is factorized according to the Brink-Axel hypothesis 
@0by 

P{E,E^) p{E ~ E^)a{E^). (1) 

The assumptions and methods behind the factorization 
of this expression are described in Refs. jl^Jl^, and only 
a short outline is given here. 

Both the level density p and the 7-energy dependent 
function cr arc unknown. In the new iteration procedure 
[ p2[ the first p° function is simply taken as a straight line, 
and the first a'^ is calculated from Eq. (1). Then new p 
and a functions are analytically calculated by minimizing 
the least square fit to the data set P. This procedure is 
repeated until a global minimum is obtained with respect 
to the values at all {E, E^) pairs. About 50 iterations 
are necessary for fitting the ^ 150 free parameters to the 
~ 1500 data points of P. Due to methodical problems 
in the first-generation procedure, we only use data with 
7-energies Ej > 1 MeV and excitation energies E > 2.5 
and 4.0 MeV in the odd-even and even-even isotopes, 
respectively p^ . 

It has been shown ||l2| that if one solution for p and 
a is known, it is possible to construct infinitely many 
solutions with the same using the substitution 

p{E ~ E^) ^ Aexp[a{E - E^)]p{E - E^) (2) 

and 

(7{E^) ^ Bcxp{aE^)a{Ej), (3) 

where A, B and a are arbitrary parameters. In the new 
product of these two functions, a factor Ai3exp(ai?) 
is left over, which is absorbed in P, since the sum 
X^B^ P{.E, E^) is undetermined. 

In the case of ^^^Dy, Fig. I demonstrates how the pa- 
rameters A and a are determined to obtain a level den- 
sity function (data points) with correct number of levels 
around the ground state (histogram). In addition, the 
parameters reproduce the level density calculated from 
the spacing of neutron resonances at the neutron bind- 
ing energy B„, see insert of Fig. 1. 

In the following, we concentrate only on the informa- 
tion given by the level density, which is assumed to be 
independent of particular 7-ray decay routes. 



III. ENTROPY FROM A SIMPLE PAIRING 
MODEL 

The level density^ p defines the partition function for 
the microcanonical ensemble, the latter being the ap- 
propriate one for statistical descriptions of isolated sys- 
tems such as finite nuclei. The partition function for the 
canonical ensemble is related to that of the microcanon- 
ical ensemble through a Laplace transform 

POO 

Z{f3)= dEpiE)exp{-PE). (4) 
Jo 

Here we have defined (3 = f/ZcsT, where T is the tem- 
perature and kB is Boltzmann's constant. Since we will 
deal with discrete energies, the Laplace transform of Eq. 
(^) takes the form 

Z(/3)-^A£;p(i?)exp (-/?£;), (5) 

E 

where AE is the energy bin used. 

In nuclear and solid state physics, thermal properties 
have mainly been studied in the canonical and grand- 
canonical ensemble. In order to obtain the level density, 
the inverse transformation 

-| />ioo 

p{E)^— dpZ{f])exp{PE), (6) 

is normally used. Compared with Eq. (^, this transfor- 
mation is rather tricky to perform since the integrand 
exp (P E + In Z ((3)) is a rapidly varying function of the 
integration parameter. In order to obtain the density of 
states, approximations like the saddle-point method, viz., 
an expansion of the exponent in the integrand to second 
order around the equilibrium point and subsequent inte- 
gration, have been used widely. For the ideal Fermi gas, 
this gives the following density of states 

cxp^2\/^ 

Pidoal(-C') - ^ r— , (7) 

where a in nuclear physics is a factor typically of the or- 
der a = A/8 with dimension MeV~^, A being the mass 
number of a given nucleus. 

Ideally, the experiment should provide the level density 
as function of excitation energy and thereby the 'full' par- 
tition function for the microcanonical ensemble. In the 



^Hereafter we use p for the level density in the microcanoni- 
cal ensemble. Furthermore, since we are dealing with a system 
with discrete energies from a quantal system, the microcanon- 
ical partition function is defined by the number of states at a 
given energy E 
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microcanonical ensemble we could then extract expecta- 
tion values for thermodynamical quantities like temper- 
ature T, or the heat capacity C. The temperature in the 
microcanonical ensemble is defined as 



(T) 



(8) 



It is a function of the excitation energy, which is the rel- 
evant variable of interest in the microcanonical ensem- 
ble. However, since the extracted level density is given 
only at discrete energies, the calculation of expectation 
values like T, involving derivatives of the partition func- 
tion, is not reliable unless a strong smoothing over en- 
ergies is performed. This case is discussed at large in 
Ref. ||] and below. Another possibility^ is to employ the 
Laplace transformation of Eq. in order to evaluate 
various thermodynamical quantities in the canonical en- 
semble. As an example, we can evaluate the entropy in 
the canonical ensemble using the definition of Helmholtz 
free energy 



F{T) = -kBT\nZ{T) = {E{T)) -TS{T). 



(9) 



Note that the temperature T is now the variable of inter- 
est and the energy E is given by the expectation value 
(E) as function of T. Similarly, the entropy S is also a 
function of T. 

In this section, we extract the exact level density from 
a simple theoretical modeP. The Hamiltonian we use to 
obtain the eigenvalues and the level density p{E) is the 
simple pairing Hamiltonian 



(10) 



ij>0 



where and a are fermion creation and annihilation op- 
erators, respectively. The indices i and j run over the 
number of levels L, and the label i stands for a time- 
reversed state. The parameter G is the strength of the 
pairing force while Si is the single-particle energy of level 
i. 

We assume that the single-particle levels are equidis- 
tant with a fixed spacing d such that Eq. becomes 



H 



2 ^ 

y>0 
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(11) 



Moreover, in our simple model, the degeneracy of the 
single-particle levels is set to 2 J -I- 1 = 2, with J = 1/2 
being the spin of the particle. 

Introducing the pair-creation operator 



and 



S 



(12) 



(13) 



one can rewrite the Hamiltonian in Eq. (|11|) as 

H^dJ2^N,-lGj2s+Sr^ (14) 



ij>0 



where 



(15) 



is the number operator. The latter commutes with the 
Hamiltonian H. In this model quantum numbers like se- 
niority S are good quantum numbers, and the eigenvalue 
problem can be rewritten in terms of blocks with good 
seniority. Loosely speaking, the seniority quantum num- 
ber S is equal to the number of unpaired particles, see 
Ref. [^g) for further details. 

The reason why we focus on such a simple model is 
twofold. Firstly, we expect the ground state of nuclei 
to be largely dominated by pairing correlations. This is 
mainly due to the strong singlet ^5*0 state in the nucleon- 
nucleon interaction, see e.g., Refs. pO|-p^. For even-even 
systems this is typically reflected in an energy gap be- 
tween the ground state and the first excited state, a gap 
that is larger than that seen in odd nuclei. This is taken 
as an evidence of strong pairing correlations in the ground 
state. More energy is needed in order to excite the sys- 
tem when all fermions are paired, i.e. when we have a 
system with an even number of particles. Since pairing 
correlations are important in nuclear systems, we expect 
the above model to exhibit at least some of the properties 
seen in finite nuclei. 

Secondly, for particle numbers up to iV ~ 20, the above 
model can be solved exactly through numerical diagonal- 
ization, since seniority is a good quantum number. This 



■^The transformation to the canonical ensemble represents 
also a smoothing. 

similar analysis within the framework of several BCS 
ansatz based approaches was done by D0ssing et al. in Ref. Q . 
Whereas our approach includes all possible eigenvalues in or- 
der to determine the level density, D0ssing et al. perform their 
diagonalization within a space spanned by number-projected 
states. The qualitative behavior of their results is however 
similar to that presented here. 
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means that we can subdivide the full eigenvalue problem 
into minor blocks with given seniority and diagonalize 
these separately. In our case we use for the even sys- 
tem = 12 particles which are distributed over L = 12 
two-fold degenerate levels giving a total of 

(^^) = (12) =2-704.156 (16) 

states. Of this total, for 5 = 0, i.e. no broken pairs, we 
have 

( N/2 ) - ( 6 ) = 924, (17) 

states. Since the Hamiltonian docs not connect states 
with different seniority S, we can diagonalize a 924 x 924 
matrix and obtain all eigenvalues with 5 = 0. Similarly, 
we can subdivide the Hamiltonian matrix into S = 2, S 
= 4,... and 5=12 (all pairs broken) blocks and obtain all 
2.704.156 eigenvalues for a system with L = 12 levels and 
TV = 12 particles. As such, we have the exact density of 
levels and can compute observables like the entropy, heat 
capacity etc. This numerically solvable model enables 
us to compute exactly the entropy in the microcanonical 
and the canonical ensembles for systems with odd and 
even numbers of particles. In addition, varying the level 
spacing d and the pairing strength G, may reveal features 
of e.g., the entropy which are similar to those of the ex- 
perimentally extracted entropy. Recall that the nuclei 
studied represent both even-even and even-odd nucleon 
systems. 

Here we study two systems in order to extract differ- 
ences between odd and even systems, namely by fixing 
the number of doubly degenerated single-particle levels to 
L = 12, whereas the number of particles is set to iV = 11 
and N = 12. Fig. || pictures the ground state of these 
systems and possible excited states. 

These two systems result in a total of ~ 3 x 10^ eigen- 
states. In the calculations, we choose a single-particle 
level spacing of d = 0.1, which is close to what is ex- 
pected for rare earth nuclei. In this sense, if we are to as- 
sign energies with dimension MeV, our results may show 
some similarity with experiment. We select three val- 
ues of the pairing strength, namely G = 1,0.2 and 0.01, 
resulting in the ratio S — d/G = 0.1, S = d/G — 0.5 
and S = d/G = 10, respectively. The first case repre- 
sents a strong pairing case, with almost degenerate single- 
particle levels. The second is an intermediate case where 
the level spacing is of the order of the pairing strength, 
while the last case results in a weak pairing case. As 
shown below, the results for the latter resemble to a cer- 
tain extent those for an ideal gas. 



*Note that the first state with a broken pair appears at a 
lower excitation energy for the odd system, as expected. 



The calculational procedure is rather straightforward. 
First we diagonalize the large Hamiltonian matrix (which 
is subdivided into seniority blocks) and obtain all eigen- 
values E for the odd and even particle case. This defines 
also the density of levels p{E), the partition function and 
the entropy in the microcanonical ensemble. Thereafter 
we can obtain the partition function Z(T) in the canon- 
ical ensemble through Eq. (|^). The partition function 
Z{T) enables us in turn to compute the entropy S(T) 
using 

S{T) = kB In Z{T) + (E{T))/T. (18) 

Since this is a model with a finite number of levels and 
particles, unless a certain smoothing is done, the micro- 
canonical entropy may vary strongly from energy to en- 
ergy. This is seen in Fig. ^, where we plot the entropy for 
the odd (upper part) and even (lower part) system using 
S = d/G — 0.5. The entropy is given by discrete points, 
since we do not have eigenvalues at all energies. How- 
ever, we can also perform a moderate smoothing which 
conserves the basic features of the model, namely an in- 
crease in entropy when pairs are broken. This was per- 
formed with a Gaussian smoothing 

- _ EkSke^pi-iE,-E,)y2a^) 
' Efeexp(-(£;,-i?fe)V2a2) ^ ^ 

where Sk and Ei^^ are the entropies and energies from 
the diagonalization of the pairing Hamiltonian. S is 
the smoothed entropy. With a smoothing parameter of 
(T = 0.2 we see that the smoothed entropy still keeps track 
of the points where the entropy experiences an increase 
due to breaking of pairs. 

Figure |^ clearly reveals the energies where two, three, 
four and so forth quasi-particles are created, i.e., where 
sudden increases in entropy take place. For the even sys- 
tem with the ground state at Eqs — —2.44, the first 
seniority 5 = 2 (formation of two quasi-particles) state 
appears at an excitation energy oi E = 2.2, the first 5 
= 4 state appears at = 4.06 and the first 5 = 6 state 
is at = 5.41. Note well that in the figures of calcula- 
tions we do not show the absolute energies. If we wish 
to employ dimensions in MeV, the first excited state for 
the even system would be close to what is expected ex- 
perimentally. 

For the odd system, the first excited states are just one- 
quasi-particle states, i.e., excitations of the last and least 
bound single-particle. Since the level spacing is much 
smaller around the ground state energy for the odd case 
(with energy Eqs = —1.65), these states appear rather 
close to the ground state. When a pair is broken, we 
create a three-quasi-particle state (one broken pair plus 
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a quasi-particle) , or seniority 5 = 3 state. This appears 
at an excitation energy'' of i? = 2.01, whilst the seniority 
5 = 5 state (two broken pairs plus one quasi-particle) 
appears at i? = 3.58. We note from Fig. || that at an 
energy of iJ ~ 8 — 9, the entropy starts decreasing (pop- 
ulation inversion), reflecting thereby the limited size of 
our model. 

For 6 = 0.5, where the single-particle spacing is only 
half the pairing strength, the energy eigenvalues are fairly 
well distributed over the given energy range. If we de- 
crease S however, we approach the degenerate limit, and 
the eigenvalues and the entropy are sharply concentrated 
around those eigenvalues where pairs are broken. This is 
seen in Fig. || for (5 = 0.1 for the even case with N = 12. 
The odd case with = 11 exhibits a similar behavior. 
Clearly, if we wish to evaluate the temperature according 
to Eq. (^) for (5 = 0.1, even with a strong smoothing, we 
cannot obtain reliable values for e.g., T. Thus, rather 
than performing a certain smoothing, we will choose to 
present further results for the entropy in the canonical 
ensemble, using the Laplace transform of Eq. (||). 

The results for the entropy in the canonical ensemble 
as functions of T for the above three sets oi 6 = d/G are 
shown in Fig. |5[ For the two cases with strong pairing, 
we see a clear difference in entropy between the odd and 
the even system. The difference in entropy between the 
odd and even systems can be easily understood from the 
fact that the lowest-lying states in the odd system in- 
volve simply the excitation of one single-particle to the 
first unoccupied single-particle state, and is interpreted 
as a single-quasi-particle state. These states are rather 
close in energy to the ground state and explain why the 
entropy for the odd system has a finite value already 
at low temperatures (recall also the discussion in connec- 
tion with Fig. H) . Higher lying excited states include also 
breaking of pairs and can be described as three-, five- and 
more-quasi-particle states. For 6 — 10, the odd and even 
systems merge together already at low temperatures, in- 
dicating that pairing correlations play a negligible role. 
For small single-particle spacing, also the difference in en- 
ergy between the first excited state and the ground state 
for the odd system is rather small. 

For our choice of d we observe that the maximum en- 
tropy is of the order of 5 ^ 14fcB in the canonical ensem- 
ble, whereas in the microcanonical ensemble, see Figs. ^ 
and ^, the maximum value is S* ^ 10 — 12kB- Obviously, 
when performing the transformation to the canonical en- 
semble, since we have a small system, there may be larger 
fluctuations in expectation values like the entropy. In the 
limit iV — > cxD, the two ensembles should result in equal 
values for T, E and S, see Rcf. ||l^ for an in depth dis- 
cussion. 



^If we wish to make contact with experiment, we could again 
assign units of MeV to ksT and E. 



For 5 = 0.5 we note that at a temperature of fc^T ~ 
0.5 — 0.6, the even and odd system approach each other'^. 
The temperature where this occurs corresponds to an ex- 
citation energy (E) in the canonical ensemble of (E) ^ 
4.7 — 5.0. Recalling Fig. ^, this corresponds to excitation 
energies where we have 4—6 quasi-particles, seniority S 
= 4 — 6, in the even system and 5—7 quasi-particles, se- 
niority iS = 5 — 7, in the odd system. The almost merging 
together of the even and odd systems at these tempera- 
tures, can be retraced to the features seen in Fig. ||. For 
higher excitation energies in Fig. ^, we saw that higher 
seniority values show less marked bumps in the entropy, 
indicating that the level density at higher excitation en- 
ergies contains many more states and that we are getting 
closer to a phase where pairing plays a less significant 
role. 

For small systems like finite nuclei, where the size of 
the system is not large compared to the range of the 
strong interaction, the entropy is not an extensive quan- 
tity, i.e., it does not scale with the size of the system 
p8[ . However, if we assume that the entropy is an ex- 
tensive quantity, then S = nSi, with n the number of 
particles and Si the single-particle entropy in the canon- 
ical ensemble. In our case Si should correspond to the 
single-quasi-particle entropy. If we label the entropy ex- 
cess AS" as the difference between the odd and even en- 
tropies, namely AS = Sodd - Seven with S'odd = «odd5'i 
and 5cvcn ~ '^cvcnS'i, we can in turn define the number 
of quasi-particles in the odd and even systems as 

nodd{E) = and novcn(£') = , (20) 

respectively. The odd system has one more quasi-particle 
than the even system, i.e., riodd = 'T-ovcn + 1- 

Fig. ^ shows the number of quasi-particles in the odd 
and even systems for the three values of S using the def- 
inition in Eq. (pO|). We note that for all cases the dif- 
ferences between the odd and even systems remain equal 
and close to one, demonstrating that the entropy is an 
extensive quantity as function of the number of quasi- 
particles. Furthermore, for S = 0.5 (central panel), we 
see that the excitation energies where 1,2,3,... quasi- 
particles appear, agree with the results discussed in Fig. 
^ in the microcanonical ensemble. To give an example, 
for the odd system, three quasi-particle appear at an en- 
ergy of (E) — 1.8, which should be compared with the 
exactly calculated one in the microcanonical ensemble 
of E = 2.01. Five quasi-particles show up at (E) — 3.4, 
which again should be compared with the result obtained 
in the microcanonical ensemble of E — 3.58. The agree- 
ment for the even case is slightly worse. For S = 0.1, the 
strong pairing case, we note that more energy is needed 
in order to create 2,4,... and 3, 5, . . . quasi-particles in 
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the even and odd systems, respectively. This agrees also 
with the microcanonical result of Fig. ^. For the weak 
pairing case S = 10, higher seniority states appear al- 
ready at low excitations energies, indicating that pairing 
plays a minor role, as expected. 

Fig. ^ carries also an interesting message. If one can 
extract the number of quasi-particles as function of exci- 
tation energies, the steepness of the curve provides use- 
ful information about the relation between the single- 
particle spacing and the pairing strength. 

In summary, varying S, allows us to extract qualitative 
informations about thcrmodynamical properties such as 
the entropy and the number of quasi-particles in even and 
odd systems. Especially, two properties are worth paying 
attention to concerning the discussion in the next section. 
Firstly, for the two cases with strong pairing ((5 = 0.1 and 
S = 0.5), Fig. II tells us that at temperatures where we 
have 4 — 6 quasi-particles in the even system and 5 — 7 
quasi-particles in the odd system, the odd and even sys- 
tem tend to merge together. This reflects the fact that 
pairing correlations tend to be less important and we ap- 
proach the non-interacting case. For the weak pairing 
case, S — 10, the odd and even systems yield similar re- 
sults at much lower temperatures. In a simple model with 
just pairing interactions, it is thus easy to see where, at 
given temperatures and excitation energies, certain de- 
grees of freedom prevail. For the experimental results 
in the next section, this may not be the case since the 
interaction between nucleons is much more complicated. 
The hope however is that pairing may dominate at low 
excitation energies and that the features seen in e.g.. Fig. 
^ are qualitatively similar to the experimental ones. 

Secondly, we can read from Fig. |6|the excitation energy 
where different numbers of quasi-particles appear. With 
a realistic value for the level spacing, a comparison with 
experiment may tell us something about the strength of 
the pairing force. 

IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 

The experimental level density p(E) at excitation en- 
ergy E is proportional to the number of levels accessible 
in 7-decay. For the present reactions the spin distribution 
is centered around (J) ~ 4.4 h with a standard deviation 
of a J ~ 2.4 h Hence, the entropy^ can be deduced 
within the microcanonical ensemble, using 

S(E) = kBlnN{E)^kBln^^, (21) 

Po 



where N is the number of levels in the energy bin at en- 
ergy E. The normalization factor pQ can be determined 
from the ground state band in the even-even nuclei, where 
we have N(E) ^ 1 within a typical experimental energy 
bin of - 0.1 MeV. 

The extracted entropies for the i6i,i62j-)y i7i,i72Yb 
nuclei are shown in Figs. and In the transformation 
from level density to entropy we use Eq. (21) with po ^ 
3 MeV^^ The entropy curves are rather linear, but with 
small oscillations or bumps superimposed. The curves 
terminate around 1 MeV below their respective neutron 
binding energies due to the experimental cut excluding 
7-rays with E^ < 1 MeV. All four curves reach S' ~ 13 
ks, which by extrapolation correspond to 5 ~ 15 ks at 
the neutron binding energy _B„ . This is the maximum en- 
tropy that a nucleus in this mass region can reach before 
neutron emission. 

The calculations for odd and even systems (see Fig. 3) 
show clear increases in the entropy at the excitation en- 
ergies where Cooper pairs are broken. This behavior is 
not very pronounced in the experimental data, probably 
due to residual couplings in real nuclei. In particular, our 
pairing model excludes collective excitations, which are 
known to contribute strongly at low excitation energy. 
For ^^^Yb in Fig. 8 one can identify bumps at 1.5 MeV 
and 2.8 MeV of excitation energy, that could be inter- 
preted as increased entropy due to the breaking of two 
and four quasi-particles, respectively. 

For the odd system the valence particle (or hole) is ex- 
pected to perform blocking, and indeed the calculations 
of Fig. 3 reveal effects of smearing out the entropy struc- 
tures as function of excitation energy. The smoother ex- 
perimental entropy curves for -'^^-'^Dy and ^^^Yb (see Figs. 
and |8|) seems also evident, in particular for the ^^^Dy 
case. 

The experimental entropy of the even-odd system fol- 
lows closely the entropy for the even-even system, but 
the even-odd system has an entropy excess. The differ- 
ence of entropy in the even-odd system compared to the 
even-even system is evaluated in Fig. |^ for I6i,i62-|-)y ^^^^ 
i7i,i72y]3^ The observed entropy excesses in the 1.5 MeV 
< E < 5.5 MeV excitation region are AS" ~ 1-8(1) ks 
and — 1-6(1) fcs for dysprosium and ytterbium, respec- 
tively. The calculations of Fig. 5 show that the entropy 
excesses are due to the additional degrees of freedom im- 
posed by a valence particle (or hole) . In the center panel 
of Fig. 5, using S = 0.5, the entropy excess is AS" ^ Skg 
at ksT — — 0.3. However, the model is based on the 
density of states formed by 1/2 spin particles. The aver- 
age total spin is J '-^ 1 giving picvci ^ Pstatc/3. Hence, the 



The experiment reveals the level density and not the state 
density. Thus, also the observed entropy reveals the num- 
ber of levels. The state density can be estimated by Pstatc ~ 
(2J + l)picvGi ~ 9.8 picVGl- 
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calculated entropy excess based on level density should 
be on the average AS* '--^ (3 — \n3)kB ~ Sfcs, which is 
close to experiment. The temperature region up to 0.3 
MeV corresponds roughly to the experimental excitation 
energy region discussed, and the model thus supports 
that excitation can be described by some sort of quasi- 
particles. 

The experimental level density can be used to deter- 
mine the canonical partition function Z{T). However, 
in the evaluation of Eq. , we have to extrapolate the 
experimental p curve to ~ 40 MeV. Here, we use the 
back-shifted level density formula of Refs. PJPq] with 



P = f 



exp[2> 



L/4[/5/V' 



(22) 



where the back-shifted energy \s U — E ~ Ei and 
the spin cut-off parameter a is defined through cr^ = 
0. 0888 A^/^y/ all. The level density parameter a and the 
back-shift parameter Ei are defined by a = G.21A°-*^ 
MeV~^ and Ei = Ci + A, respectively, where the cor- 
rection factor is given by Ci = —6.6A~^-^^ according to 
Rcf. The factor / is introduced by us to adjust the 
theoretical level density to experiment at ~ B„ — 1 
MeV. The parameters employed are listed in Table 1. 
From our semi-experimental partition function, the en- 
tropy can be determined from Eq. (18). The results are 
shown in Fig. The entropy curves show a splitting at 
temperatures below fcsT — 0.5 — 0.6 MeV which reflects 
the experimental splitting shown in the microcanonical 
plots of Figs. 7 and 8. However, the strong averaging pro- 
duced by the summing in Eq. (5), modifies the entropy 
excess due to components from the theoretical extrapola- 
tion of p. Even so, the curves agree qualitatively with the 
calculations in Fig. 5 using S = 0.5. The effect of pair- 
ing seems in both cases to vanish above 0.5 — 0.6 MeV. 
This agrees with our previous work ]To| , giving a critical 
temperature of fesTc = 0.5 MeV for the existence of pair 
correlations. 

The observation that one quasi-particle carries 1.7 ks 
of entropy, can be utilized to estimate the number of 
quasi-particles as function of excitation energy. Analo- 
gously to Eq. (20), we estimate from the experimental 
entropies Sco and Sac in neighboring even-odd and even- 
even isotopes the entropy excess AS" = Seo — See- The 
number of quasi-particles in the even-odd and even-even 
systems is given by Eq. (|20|), except that the odd system 
is replaced by an odd-even nucleus and the even system 



by an even-even nucleus. 

The extracted number of quasi-particles n{E) in ^^^Dy 
and ^^^Yb is shown in Fig. |ll|. The number of quasi- 
particles raises to a level of n ^ 2 around E — 1.5 — 2 
MeV, which could be a signal for the formation of two 
quasi-particle states. However, the creation of four and 
six quasi-particles shows no clear step-like function. The 
breaking of additional pairs is spread out in excitation 
energy giving a rather smooth increase in the number 
of quasi-particles as function of excitation energy. In the 
excitation region 0.5 — 5 MeV the n{E) curve gives on the 
average 1.6 MeV of excitation energy to create a quasi- 
particle pair. This value is consistent with pairing gap 
parameters of this mass region, see Table 1. The theo- 
retical calculation'' with S = 0.5 gives an energy of 1.7 
MeV per broken pair, which is close to the experimental 
finding of 1.6 MeV. Hence, with a single-particle spacing 
of d = 0.1 — 0.2 MeV, the pairing strength is determined 
to G = 0.2 - 0.4 MeV. 



V. CONCLUSIONS 

The entropy as function of excitation energy has been 
extracted for the ^^^''^^'^Dj and i7i,i72Yb nuclei. The ob- 
served entropy excess in the even-odd nuclei compared 
to the even-even nuclei is interpreted as the entropy for 
a single quasi-particle (particle or hole) outside an even- 
even core. The entropy excess remains at a level of AS* ~ 
1.7 kB as function of excitation energy. A simple pairing 
model with an equidistant level spacing of d and a pairing 
strength of G, gives a qualitatively similar description of 
these features. 

The number of excited quasi-particles has been ex- 
tracted from data. The onset of two quasi-particle ex- 
citations seems evident; however, the breaking of addi- 
tional pairs is smeared out in excitation energy and is 
difficult to observe. The maximum number of excited 
quasi-particles is measured to n ~ 6 at an excitation en- 
ergy of 5.5 MeV in the ^^^Dy and ^''^Yb isotopes. 

The quasi-particle picture has been a success in de- 
scribing rotational bands in cold nuclei. The present re- 
sults indicate that quasi-particles also can describe cer- 
tain thermodynamical properties of hot nuclei. This gives 
hope for realistic modeling of nuclei up to high intrinsic 
energy with several quasi-particles excited. 

The authors are grateful to E.A. Olsen and J. Wikne 
for providing the excellent experimental conditions. We 



''The reader should keep in mind that the number of parti- 
cles in the theoretical calculation and experiment are rather 
different. In experiment, if one assumes ^'''^Sn as closed shell 
core, the number of valence protons and neutrons is of the or- 
der of ~ 30 — 40. However, performing the above theoretical 
calculations with say 10 or 14 particles results in qualitatively 
similar results as those presented here. 
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Table 1: Parameters used in the back-shifted Fermi gas formula for the extrapolation of the experimental level density 
curve. 



Nucleus 


A (keV) 


a (MeV-^) 


Ci (keV) 


El (kcV) 


/ 


i«Dy 


793 


17.46 


-1298 


-505 


1.400 




1847 


17.56 


-1296 


551 


1.138 


17lYb 


680 


18.40 


-1273 


-593 


0.376 


172Yb 


1606 


18.50 


-1271 


335 


0.465 
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FIG. 1. The extracted level density for Dy reproduces both known levels (histogram) at low excitation energy and the 
neutron resonance spacing at B„ (triangle) by proper adjustment of the parameters A and a of Eq. (2). 
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FIG. 2. Simple illustration of the ground state and possible excited states ( 4 and 5 quasi-particles) for a system with 

12 doubly-dcgcnorato singlc-particlo levels. The properties of the model is governed by the level spacing d and the pairing 
strength parameter G (the illustration is with G = 0). For the even system with 12 particles, the first excited state is a two 
quasi-particle state corresponding to the breaking of one pair. The first excited state in the odd system with 11 particles is a 
single quasi-particle state. 
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FIG. 3. Entropy in the microcanonical ensemble as function of excitation energy E for 5 = 0.5. The upper and lower panels 
show the results for the odd and even systems, respectively. Results with and without a Gaussian smoothing are displayed. 
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FIG. 4. Entropy in the microcanonical ensemble as function of excitation energy E for S = 0.1. 
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FIG. 5. Entropy in the canonical ensemble as function of temperature ksT for odd and even systems for S = 0.1 (upper 
panel), 6 = 0.5 (central panel) and 5 = 10 (lower panel). 
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FIG. 6. Number of quasi-particles n in the canonical ensemble for different values of S for even and odd particle systems. 
Results for d = 0.1 are shown in the upper panel, S = 0.5 in the central panel and 5 = 10 in the lower panel. 
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FIG. 7. Observed entropy for ' Dy as function of excitation energy E. 
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FIG. 8. Observed entropy for ' Yb as function of excitation energy E. 
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FIG. 9. The entropy excess /S.S in Dy compared to Dy (upper panel) and in Yb compared to Yb (lower panel). 
The lines through the data points indicate the average values found. 
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FIG. 10. The semi-experimental entropy S for ' Dy and ' Yb calculated in the canonical ensemble as function of 
temperature fcsT. 



18 



c 

O 

o 



1 62 



4. ^ * f .1* * 



X3 



3 



_C1J 
O 



o 

CL 




12 3 4 

Excitation energy E [MeV] 



FIG. 11. The number of quasi-particles n in Dy (upper panel) and Yb (lower panel) as function of excitation energy. 
The lines indicate the levels of two, four and six quasi-particles. 
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